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ANNOTATED  COMPUTER  OUPUT  FOR  SPLIT  PLOT  DESIGN:  BMDP  2V 

by 

W.T.  Federer,  Z.D.  Feng,  and  N.J.  Miles -McDermott 
Mathematical  Sciences  Institute,  Cornell  University,  Ithaca,  NY 

ABSTRACT 

In  order  to  provide  an  understanding  of  covariance  analyses  for 
split  plot  designs,  analyses  are  first  conducted  on  a  data  set 
without  the  covariate.  Then  an  analysis  with  the  covariate  is 
performed.  The  third  example  has  a  different  experiment  design  for 
the  whole  plots  and  the  covariate  is  constant  for  all  split  plots 
within  a  whole  plot.  Computer  packages  may  only  give  portions  of 
the  analysis  correctly.  Some  packages  require  two  procedural  calls 
to  obtain  a  portion  of  the  correct  results.  BMDP  2V  falls  in  the 
above  category. 

INTRODUCTION 

This  is  part  of  a  continuing  project  that  produces  annotated 
computer  output  for  the  analysis  of  balanced  split  plot  experiments 
with  covariates.  The  complete  project  will  involve  processing  three 
examples  on  SAS/GLM,  BMDP/2V,  SPSS-X/MANOVA,  GENSTAT/ANOVA ,  and 
SYSTAT/MGLH.  Only  univariate  results  are  considered.  We  show  here 
the  results  from  BMDP  2V. 

For  Example  1,  the  data  are  artificial  and  were  constructed 
for  ease  of  computation;  the  experiment  design  for  the  whole  plots 
is  a  randomized  complete  block  and  the  split  plot  treatments  are 
randomly  allocated  to  the  split  plot  experimental  units  within  each 
whole  plot.  Example  2  is  the  same  as  Example  1  except  that  a 


covariate  varies  from  split  plot  to  split  plot.  The  data  for 
Example  3  come  from  an  experiment  wherein  the  whole  plot  treatments 
were  laid  out  in  a  completely  randomized  design  and  the  split  plot 
treatments  were  randomly  allotted  to  the  split  plot  experimental 
units  within  each  whole  plot.  The  value  of  the  covariate  varies 
from  whole  plot  to  whole  plot  but  is  constant  for  all  split  plots 
within  a  whole  plot  treatment. 

We  present  the  elementary  computational  steps.  Simple 
hypothetical  data  are  used  for  the  first  two  examples  so  that  it  is 
easy  to  provide  all  detailed  computations  to  illustrate  how  each 
number  is  obtained.  Some  readers  may  wish  to  skip  the  detailed 
computations.  The  third  example  comes  from  Winer  (1971) .  The 
detailed  computations  are  given  in  his  book  (p.  803) . 

Data  SP-1 

Split  plot  data  with  whole  plots  arranged  in 
randomized  complete  block  design 
(hypothetical  data) 


W1 

Whole  plot 

treatment 

W2 

split 

S1 

plot 

S2 

treatment 
S3  S4 

Total 

split 

S1 

plot 

S2 

treatment 
S3  S4 

Total 

3 

4 

7 

6 

20 

3 

2 

1 

14 

20 

6 

10 

1 

11 

28 

8 

8 

2 

18 

36 

6 

10 

4 

4 

24 

10 

8 

9 

13 

40 

15 

24 

12 

21 

72 

21 

18 

12 

45 

96 

Block 


I 


Total  and  Means 


Blocks  W(whole  plots)  S (split  plot) 

(8  observations)  (12  observations)  (6  observations) 


Mean 

Mean 

Mean. 

1 

40 

5 

W1  72 

6 

Sx  36 

6 

2 

64 

8 

W2  96 

8 

S2  42 

7 

3 

64 

8 

S3  24 

4 

Grand 

Total 

168 

S4  66 

11 

Grand 

Mean 

7 

Model:  *ijk 

■  4  +  Pj  +  + 

6ij +  ak 

+  (OT)lk  +  Cijk 

4 

mean 

T  i 

=  effect 

of  whole  plot  i 

Pj  = 

effect  of 

block  j 

=  effect 

of  split  plot  k 

6ij  = 

error  (a) 

effect  <aT ) iR 

=  effect 

of  interaction  of 

€ijk 

error  (b) 

effect 

whole  plot  i  and  split 

plot 

Analysis  of  Variance 


Source 

(*) 

df 

SS 

B  (Blocks) 

« 

R(p 

4  , r ,a ,ar ) 

2 

48 

W  (whole  plot  treatments) 

= 

R(t 

4 , p ,a ,ar ) 

1 

24 

BxW  (error  (a)) 

= 

R(6 

4,p, T,a,aT) 

2 

16 

S  (split  plot  treatments) 

= 

R(a 

4,p ,T ,ot) 

3 

156 

SxW  (interaction  of  S  and  W) 

ss 

R  (ar |4,a,r,p) 

3 

84 

(**)  SxB: W  Terror  fbn 

= 

R(€ 

li.a.T.aT.p) 

12 

Total  (Corrected  for  mean) 

= 

R(p 

r,6,a,aT,€  p) 

23 

440 

Mean 

= 

R(4) 

1 

1176 

Total  TUncorrected  for  meanl 

= 

.  Rf-^-j 

.Pi ti6x°i°T,€) _ 

_ ZA _ 

(*)Notation  follows  that  of  Searle ( 1971) .  Since  the  design  is  balanced, 
R(P  |4  ,a  ,ar)  =  R(p|fi),  etc.  The  simpler  notation  is  used  later. 

(**)  SxB:W  means  SxB  within  W. 


Calculations  of  sums  of  squares: 

N  *  2*3*4  =24  ,  Y  =  7 

2  3  4 

R(4,p,T,6,a,aT,€)=  222  Y?.v  =  (32+  62+  62+***+  182+  132)  =  1616 

i=l  j=l  k=l  1]K 

R(4)  =  NY2  =  24 • (7 ) 2  *  1176 
R(p,T,6,a,aT,e\»)  =  1616  -  1176  =  440 


3 


2  2  2 

R(p|n)  -  R (*i,P)  -  R(m)  =  (4°  ±*±J±$±J-  -  H76  =  1224  -  1176  =  48 

2  2 

R(t|m)  =  ROi,T)  -  R  (p )  =  1296  ^  ”  1176  *  1200  -  1176  =  24 

R(6|m»P#t)  =  R(6  ,n  ,  p  ,t  )  -  R(p,p)  -  R(t,»i)  +  R  (ji ) 

=  .  1224  .  1200  +  11?6 
4 

=  1264  -  1224  -  1200  +  1176  =  16 


R(a  | m )  =  R (a,\i)  -  R(fi) 


2  2  2  2 
=  (36Z+42‘*+24^+66<g) 

6 


1176  =  1332  -  1176  =  156 


R(aT|n,a,T)  *  R(ar,4,a,T)  -  R(fi  ,a)  -  R(p,r)  +  R(p) 

-  (lS2+2«2+122+2lW+182+12W)  .  ,33,  .  1200  +  11?6 

=  1440  -  1332  -  1200  +  1176  =  84 

R(t  Im  ,P  ,6  ,a  ,t  ,ar)  *  R  (t  ,u  ,a  ,  p  ,  6  ,r  ,ar )  -  R(p,p,T,6)-  R(p,a,T,aT)  +  R(t,ji) 

*  1616  -  1264  -  1440  +  1200  *  112 


Data  SP-2 

Data  SP-2:  Data  SP-1  with  the  following  covariate  Z 
which  varies  with  split  plot 


Covariate  (Z) 


whole  Dlot 

W1 

L  W2 

S1 

S2 

S3 

S4 

Total 

S1 

S2 

S3 

S4 

Total 

B1 

1 

2 

1 

2 

6 

2 

0 

2 

■ 

8 

B2 

2 

2 

0 

4 

8 

4 

1 

3 

12 

B3 

3 

5 

2 

0 

10 

3 

2 

4 

9 

16 

Total 

6 

9 

3 

6 

24 

9 

3 

9 

15 

36 

Totals  and  Means 


blocks 

(8  observations) 
otal  Mean 


1 

14 

14 

2 

20 

20 

3 

26 

26 

W  (whole  plot) 
(12  observations) 
otal  Mean 


S  (split  plot) 
(6  observations) 
Total 


l 

l 

1 

2 


Model:  Yijk  -  M+pj+Ti+6i;j+ek+(«T).k+Pl(2ij.-Z...)+P2(Zijk-2ij.)+tijk 


pj  =  effect  of  jth  block 

=  effect  of  kth  split  plot 
/3  =  split  plot  regression  slope 


=  effect  of  ith  whole  plot 
P1  =  whole  plot  regression  slope 


=  error  a 


‘ijk  -  error  b 


Table  of  sum  of  squares  and  products 


Source _ 

B 

W 

BxW  (error  a) 

S 

SxW 

SxB:W  (error  b) 

Mean _ 

Total _ 


YY 

48 

24 

16 

156 

84 

112 

l!76 

.616 


YY  column  is  the  same  as  in  SP-1,  ZZ  column  is  computed  in  the  same 
fashion.  Thus,  only  computations  for  YZ  column  are  illustrated. 


Total  =  2  2  2  Y..  -Z,. 

YZ  i=l  j=l  k=l  1DK 

=  3(1)+  6(2)+***+  14(4)+  18(4)+  13(7) 

Meanyz  =  NY. . .  Z. . .  =  «  420 


2  4 

(  2  2 
i=l  k=l 


Yijk 


2  4 

)(  2  2  Z..k) 

i=l  k=l  13K 


40(14) +64 (20)+64 (26 


4 

2 


SiRv 

W 


432  +  420 


MS 

!« 


432  +  420 


Mm 

0 

m 

$ 


2  2(2  Y..  ) (  2  Z..  ) 

i=l  k=l  t=1  t=1  1 J * 

SxWv„:  —  ■— — -  453  -  432  +  420 

YZ  3 

=  498  -  453  -  432  +  420  *  33 


SxB:Wyz: 


2  3  4 

2  2  2  Y.  ..Z.  ..  -  454  -  498  +  432 

1=1  j-l  k=l  13k  i]k 

=  537  -  454  -  498  +  432  =  17 


Source 


B  (block) 

W  (whole  plot  treatment) 
Regression  (a) 

BxW  (error  (a)) 

S  (split  plot  treatment) 

SxW  (interaction  of  S  and  W) 
Regression  (b) 

SxB:  W  (error  (b) ) 


R(P ,t) 

R(t  |m,p,/31) 

R(P1U/P#t) 

R(6  |pi  ,p  ,r  ,PX) 

R(a\n,p  ,t  ,ar  ,P2) 

R(ot \n,p ,t ,a,P2) 
R(P2\u,p ,t ,a,aT) 

R(t  \u , p ,a ,t ,qt ,P~) 


3.4286 

16.0 

0.0 

84.243 

37.474 

14.450 

97.550 


m 


Px  -  BxWyz  /  BxWzz  =4/1=4 


P2  =  SxB:Wyz  /  SxB:Wzz  =  17/20  =  0.85 


The  SS's  adjusted  by  regression  on  Z  are  illustrated  below: 


R(p|p)  =  48,  remains  same  since  it  is  not  of  interest  to  adjust  for  Z  on 
the  blocks. 


R(t,6  =  (Wyy+  BxWyy)  - 


(WYZ+  BxWyz> 
wzz+  BxWzz 


256 


-  (24  +  16)  -  H2  +  4>  =  40  -  =  3.4286 

6  +  1  7 


R(6  |  U,P,T  ,PX)  *  BxW 


YY 


-  16  -  li  .  0 


BxW 


ZZ 


Rf+lu^P#^)  =  R(t,6  |fi,p  ,P1)  -  R(6  |m/P#t,P1) 


=  40  -  -  0  =  3.4286 


■  (B*Wyz>  4 2 

R(P1lfx,T,p)  =  - =  —  =  16 


BxW 


ZZ 


R(«.e \v,p  ,t  ,clt  ,P0)  =  (Syy  +  SxB:Wvv)  - 


YY; 


(Syz  +  SxB:Wyz) 

szz  +  SxB:Wzz 


=  (156  +  112)  -  12211221 

9+20 

=  268  -  86.207  =  181.793 


(SxW  +  SxB:W  ) 

R(ar,t  |M,p,a,T,/35)  =  (SxWyy  +  SxB:Wyy) - 


SxWzz  +  SxB:WZZ 


=  84  +  112  -  =  196  -  60.976  =  135.024 

21+20 


Note:  R(a,e  |fx,/32)  and  R(oct  ,e  |p  ,a  ,t  ,/?2)  are  intermediate  steps  for 

later  use. 


7 


tl 

8# 

?$ 


«SV, 

»»Iv, 

itiV 

Vi  V 
W 

1 


£ 

•JKJ 

i| 


R02ln,p,a,T,aT)  = 


(SxB:Wyzr  _  1?2 

SxB:Wzz  20 


=  14.450 


(SxB:W  )*  172 

R(e|fi,p,a,T,aT,/3  )  =  SxB:W - — -  =112  -  =112  -  14.45  =  97.55 

XX  SxB:Wzz  20 

R(a  |p  ,  p  ,t  ,<xt  ,P2)  =  R(“  f £  ||i  »P  »t /aT  i  $2^  ~  error  b  =  181.793  -  97.55 

=  84.243 

R(ar  \u,p,a,T,P2)  =  R  (aT  ,  e  |  v ,  p  ,  a  ,  t  ,  P2 )  -  R(e  |p  ,  p  ,a  ,  t  ,«t  ,P2) 

=  135.024  -  97.55  =  37.474 
Data  SP-3 

Split  plot  data  with  plots  arranged  in  a  completely 
randomized  design  and  a  covariate  Z  that  is  constant 
within  the  whole  plot.  (Winer,  1971,  p.  803) 


:  y-  •,  =  u  +  t.  +  6..  +  a,  +  (ra)  ..  +  8 , ( Z . . -Z  )  +  e... 

Jfi}k  K  l  13  k  v  7ak  Kl'  13  13k 


=  A  effect  (whole  plot) 
=  B  effect  (split  plot) 


5 i j  =  error  (a)  e ^  =  error  (b) 
P1  =  whole  plot  regression  slope 


Analysis  of  variance  and  covariance 


Source _ _ df _ SS 


A  (whole  plot) 

=  r(t  |p.  ,/31) 

1 

44.492 

Regression 

=  R(/31|m,t) 

1 

166.577 

Error  (a) 

=  R(6  ||* ,T ,PX) 

5 

61.298 

B  (split  plot) 

=  R(a|n,T,crr) 

1 

85.563 

AxB  (interaction) 

=  R(ra  |jx  ,t  ,a) 

1 

0.563 

Error  (b) _ 

=  R(& lu ,t ,a,Ta) _ 

6 

6.375 

Table  of  SS 

and  products 

Symbol 

Y2 

ZY 

Z2 

W 

68.06 

12.38 

2.25 

E  (a) 

227.88 

163.00 

159.50 

S 

85.563 

0 

0 

WS 

0.563 

0 

0 

E(b) _ 

6.375 

0 

0 

A 


163.00 

159.50 


1.02 


Since  the  computations  are  illustrated  in  Winer  (1971,  p.  803-5)  we  have 
omitted  them  here. 


SP-1:  Control  Language 

Control  Language  is  typed  upper  case  and  comments  are  in  bold  and  lower 
case. 

/PROBLEM  TITLE  IS  'SP-l:  WHOLE  PLOT  ARRANGED  IN  RCB' . 

/INPUT  VARIABLES  ARE  6. 

FORMAT  IS  FREE. 

/VARIABLE  NAMES  ARE  BLOCK,  WHOLE,  Yl,  Y2,  Y3,  Y4 .  =*  input  variables 
/DESIGN  GROUPING  ARE  BLOCK,  WHOLE.  ==>  whole  plot  classification 

variables 

DEPENDENT  ARE  Yl,  Y2,  Y3,  Y4 .  =>  specifies  Y  measured  for 

four  levels  of  split  plot 

LEVEL  IS  4.  =>  specifies  four  levels  of  split  plot 

nAME  IS  SUBPLOT.  =>  split  plot  classification  variable 
EXCLUDE  IS  12.  =$  model  specification  for  whole  plot:  assumes 

full  factorial.  This  statement  requests  the 
BLOCK  (1)  by  WHOLE  (2)  interaction  to  be 
pooled  with  SSE(a) 

PRINT.  =>  requests  residuals  and  predicted  values 

/END 

1  1  3  4  7  6  data  organized  with  all  split  values  for  a  partic- 

2  1  6  10  1  11  ular  BLOCK  by  WHOLE  combination  on  the  same  line 

3  1  6  10  4  4 


SP-2:  Control  Language 

Note:  The  whole  plot  and  split  plot  analysis  must  be  requested  in 

separate  runs.  The  first  run  gives  the  correct  whole  plot  analysis  but  an 
incorrect  split  analysis  (details  follow  in  the  annotations) .  The  second 
run  gives  the  correct  split  plot  analysis. 

Whole  Plot  Run 


/PROBLEM  TITLE  IS  'SP-2:  SP-1  WITH  COVARIATE  VARYING  WITH  SPLIT  PLOT'. 
/INPUT  VARIABLES  ARE  10. 

FORMAT  IS  FREE. 

/VARIABLE  NAMES  ARE  BLOCK,  WHOLE,  Zl,  Z2,  Z3,  Z4,  Yl,  Y2 ,  Y3,  Y4 . 
/DESIGN  GROUPING  ARE  BLOCK,  WHOLE. 

DEPENDENT  ARE  Yl,  Y2 ,  Y3,  Y4 . 

COVARIATE  ARE  Zl,  Z2,  Z3,  Z4.  =>  yields  an  analysis  of 
LEVEL  IS  4 .  covariance 

NAME  IS  SUBPLOT. 

EXCLUDE  IS  12. 

RESIDUAL  IS  MEAN. 

PRINT. 

/END 

1112123476 
2122046  10  1  11 
3135206  10  44 
122024321  14 
224134882  18 
323247  10  89  13 


Split  Plot  Run 


/PROBLEM 

/INPUT 

/VARIABLE 

/DESIGN 


/END 

11113 

112  2  4 

113  17 

114  2  6 


TITLE  IS  'SP-2:  SP-1  WITH  COVARIATE  VARYING  WITH  SPLIT  PLOT'. 
VARIABLES  ARE  5. 

FORMAT  IS  FREE. 

NAMES  ARE  BLOCK,  WHOLE,  SUBPLOT,  Z,  Y. 

GROUPING  ARE  BLOCK,  WHOLE,  SUBPLOT. 

DEPENDENT  IS  Y. 

COVARIATE  IS  Z. 

EXCLUDE  IS  13,123.  =>  model  specification  for  split  plot: 

PRINT  program  assumes  full  factorial.  This 

statement  specifies  certain  factors 
to  be  pooled  in  SSE(b). 


for  split  plot  analysis  data  must  be  organized  in  a 
similar  arrangement  with  only  one  datum  per  line 
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SP-3 :  Control  Language 


Note:  This  example  unlike  the  other  two  covariate  examples  requires  only 
one  run.  This  is  because  we  are  using  a  CRD  and  no  special  design 
specification  is  required. 


/PROBLEM 

/INPUT 

/VARIABLE 

/DESIGN 


/END 
1  3  10  8 
1  5  15  12 

1  8  20  14 
1  2  12  6 

2  1  15  10 
2  8  25  20 
2  10  20  15 
2  2  15  10 


TITLE  IS  ' SP-4 :  WHOLE  PLOTS  ARRANGED  IN  CRD  WITH  CONSTRAINT 

COVARIATE 

VARIABLES  ARE  4. 

FORMAT  IS  FREE. 

NAMES  ARE  A,  Z,  Yl,  Y2 . 

GROUPING  IS  A. 

DEPENDENT  ARE  Yl,  Y2 . 

COVARIATE  IS  Z,  Z.  =$  BMDP  requires  that  two  covariables  be 
LEVEL  IS  2 .  specified  here  to  match  the  number  of 

NAME  IS  B.  levels  of  the  split  plot.  Since  Z  is 

PRINT.  constant  for  both  levels,  its  name  is 

entered  for  both. 
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